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Abstract 

We analyse second order (in Riemann curvature) geometric flows (un-normalised) on locally ho- 
mogeneous three manifolds and look for specific features through the solutions (analytic whereever 
possible, otherwise numerical) of the evolution equations. Several novelties appear in the context 
of scale factor evolution, fixed curves, phase portraits, approach to singular metrics, isotropisation 
and curvature scalar evolution. The distinguishing features linked to the presence of the second 
order term in the fiow equation are pointed out. Throughout the article, we compare the results 
obtained, with the corresponding results for un-normalized Ricci fiows. 
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I. INTRODUCTION AND OVERVIEW 

Ricci and other geometric flows were introduced in mathematics by Hamilton 

and in physics, by Friedan around almost the same time, though with very different 
motivations. More recently, such geometric flows have become popular, largely because 
of Perelman's work which led to the proof of the well-known Poincare conjecture. In 
physics, geometric flows have been investigated in varied contexts such as general relativity, 
black hole entropy {gI and string theory E^SL 



The unnormalised Ricci flow equation 



dgij 



2| is given as, 

= -^R^J (1) 



dt 

where Qij is the metric tensor, Rij the Ricci tensor and t denotes a time parameter (not the 
physical time). For a normalised Ricci flow, the corresponding equation turns out to be: 

^ = -2R., + l{R)g., (2) 

where (R) = "j^r- Thus the normalised flow ceases to be different from the unnormalised 
one if we consider non-compact, infinite volume manifolds (with a finite value for J RdV). In 
addition, for constant curvature manifolds, the second term in the normalised flow equation 
reduces to \Rgij- 

Friedan's early work showed how one may arrive at the vacuum Einstein field equations 
of General Relativity from the Renormalisation Group (RG) analysis of a nonlinear a-model 



a 



1^. The RG /3-function equations (after a a' dependent scaling of the flow parameter) are 
given (upto third order in the a' parameter) as follows 

dg 



dt 



-2Rc - a'Rc - 2a'^Rc - ... (3) 



where Rc is a symmetric 2-tensor defined as: Rc = RikimRjabcg^"' g''^ g"^'^ and Rc = 
\RkimpRi^^'' Rj^^ ^ - \RikijR^'''^'^ R'' spr + In component form, we have, 

= —2Rij — a RikimRj ^ — . . . (4) 

which may be considered as a higher order geometric flow equation for the metric on a given 



manifold. We intend to continue our investigations ll| on such higher (mainly second) 



order geometric flow equations, in this article. We mention that apart from a motivational 



perspective, our work has no connection with RG flows and, in what follows, we will not 
discuss any link with them anywhere, a' is treated as a parameter without any link to the 
RG flow equations. In other words, the higher order flows we talk about are treated as 
geometric flows in their own right. Restricting to second order we mention a scaling feature 
of this geometric flow. For Ricci flow a scaling of the metric Qij ^Qiji can be compensated 
by a scaling of t by t' = ^. Since the Ricci is invariant under such scaling, flows for Qij and t 
are equivalent with those for g[- and t' . For the second order flows, things are a bit different. 
Here, we can see that flows for Qij, t and a' will be equivalent to flows for Qij/Q, t' = t/Q 
and a" = a'/Q. In our work, we have retained the a' when we do analytic calculations 
-however, for numerics we have chosen values of a' = 0, ±1- solutions for other values of a' 
can be obtained using the scaling property mentioned above. 

In this paper, we investigate higher (second) order flows on locally homogeneous three 
manifolds. After the seminal work of Isenberg and Jackson [l^ on the behavior of Ricci flow 
on homogeneous three manifolds, several interesting papers on this topic have appeared in 



the literature. Notable among them are the articles in 1^, [l^, [isl, jlo ]. 



In particular. 



the authors in 13|] discuss a new approach wherein the Ricci flow of left-invariant metrics on 
homogeneous three manifolds is equivalently viewed as a flow of the structure constants of 
the metric Lie algebra w.r.t. an evolving orthonormal frame. Consequences for some other 



18|,(l9|. 



geometric flows on such three manifolds are available in 

Homogeneous spaces are also known to play an important role in cosmology in the context 
of the so-called Bianchi models ( 20[, 21[). It is worth mentioning that the Misner- Wheeler 
minisuperspace deals with three dimensional homogeneous spaces in cosmology where dif- 
ferent paths in superspace corres pon d to different evolution proflles for metrics. For more 



details in this context see 



12|, [20|, |22|, [23|. Recently, a connection between self-dual 



gravitational instantons and geometric flows of all Bianchi types has been shown in 24 1. 

Our article is organised as follows. In Section II, we discuss some preliminaries. Section 
HI reviews the general theory of left invariant metrics. Section IV contains the main results 
of the article-here we analyse the higher order flows for the various Bianchi classes. Finally, 
in Section V we conclude with some remarks. 
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II. PRELIMINARIES ON HOMOGENEOUS THREE MANIFOLDS 



Following Isenberg and Jackson 12| (to which we refer for details concerning the follow- 
ing discussion), we take the view point that our original interest is in closed Riemannian 
3- manifolds that are locally homogeneous. By a result of Singer 25|, the universal cover of 
a locally homogeneous manifold is homogeneous, i.e. its isometry group acts transitively. 



Here we begin from the basic definition of homogeneous manifolds and give Thurston's [26 1 
proposition for eight geometric structures. A homogeneous Riemannian space is a Rieman- 
nian manifold {M,g) whose isometry group lsom{M,g) acts transitively on M, i.e given 
x,y ^ M there exists a G Isom(M, g) such that = y. 

If G is a connected Lie group and H is a. closed subgroup, G/H is the space of cosets 
{gH}, TT : G — 7- G/H is defined by — )■ [gH]. G/H is called a homogeneous space. 

Let {M^, G, Gx) be a maximal model geometry i.e. the isotropy group G^ {g G G\gx = x) 
is maximal among all subgroups of the diffeomorphism group of that have compact 
isotropy groups. Depending on the isotropy group, there are eight possible manifestations 
of the maximal model geometry which are also known as Thurston geometries. Here we 
mention them briefly. If the isotropy group is 5*0(3) then the model geometry is any one 
of the three, namely, SU(2), M.^ or M.^. On the other hand, if 5*0(2) is the isotropy group, 
then it has four possibilities depending on whether can be a trivial bundle over a two- 
dimensional maximal model or not. The first possibility where is a trivial bundle over a 
2-dimensional maximal model, gives rise to x M or x M upto diffeomorphism and the 
latter produces Nil or 5*1/(2, R) up to isomorphism. The last case where the isotropy group 
is trivial, the Lie group is isomorphic to Sol. 



III. GENERAL THEORY FOR LEFT INVARIANTS METRICS ON 3D UNIMOD- 
ULAR LIE GROUPS 

In studying curvature properties of left invariant metrics on a Lie group, many results 



lave been obtained in the past. Most of them are contained in a survey article by Ml 



nor 



30|. 



27I 1. For more details on Lie groups and homogeneous spaces we refer to [28|], [29| and 
We review this work here largely because, later, we will use the results quoted here while 
writing down the second order geometric flow equations. 
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The left invariant metric on is taken as: 

g{t) = A(t)r/^ ® + B{t)r]'^ ®rf + C{t)r]^ ® r]^ 



(5) 

= tt' (g) vr' + vr^ (g) TT^ + vr'^ (g) vr'^ 

and its inverse as, 

= ei (g ei + 62 ® 62 + 63 (g) 63 

where -B(^), C{t) are positive. Here {Fj}^^]^ are the left invariant frame field (also 

called Milnor's frame) with dual coframe field {ri^}f^i . The Lie brackets w.r.t the left 
invariant frames are of this form: 

[Fi,Fj]=aijkFk (7) 

Let us define an orthonormal frame field {6j}f^i which will be of the form: 61 := "^-^i- If 
we let Ci := A, C2 := B, (3 := C and A := 0231, /i := ^312, v := 0123 then Eq.(l7]) looks like 

[ei, ej] = r—r— ek (8) 
The relevant components of the sectional curvature turn out to be 

K(6i A 62) := {Rm{eu 62)62, 6,) = ^ ^^^^ ' + ^ (9) 

K[62 A 63) := {Rm[62, 63)63, 62) = — . . + 7^7^ (10) 



AABC 


+ 


AAB 




ifiB - pCf 


^ A(2z/C- 




-3AA) 


AABC 




ABC 




{uC - XAf 


fi{2XA- 


^2vC 


-3;u5) 



K{6; A 61) := (i?m(63, 61)61, 63) = + — ^ (11) 

Recall that these are actually R12121 R3131, respectively, in the orthonormal frame. 

We have to write them back in Milnor's frame. Thus the Riemann tensor in the Milnor 
frame will be 

Rm{Fi, Fj, Fi, Fj) = (iQRm{6i, 6j, 6,, 6j) (12) 

From now on we would mean:(i?m(6j, 6j)6j, 6^) = Rm{6i,6j,6i,6j). Similarly the Ricci and 
Re would look as : 

3 1 

Rc{F„ Fi) = V Rm{F„ F^, F,, F,)- (13) 
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k=l 



2li 

k) >2 /- 



(14) 



Rc(F2,F2) 
Rc(F3,F3) 







uCf 




2BC 






- (z/C- 


\Af 




2CA 




{uCf 


- {\A- 


fiBf 



2AB 



(15) 
(16) 
(17) 



RcCF\>i) = 


1 

8AB^ 


'{XA 


-y.Bf 
C 


+ u{2^B + 2XA - 


'iuC) 


+ 


1 


'{yC 


-\Af 


+ /i(2AA + 2uC - 
+ \{2uC + 2/i5 - 


3fiB) 
3XA) 


RcCFr>2) = 


8AC^ 
1 

8BC^ 


-{liB 


B 

-vCf 
A 


+ 


1 


'[\A 




+ u{2^B + 2AA - 
+ li{2\A + 2uC - 


3uC) 
SfiB) 




1 

8A^C 


'{yC 


C 

-\AY 
B 


1 


1 


UliB 


-uCf 


+ A(2z/C + 2/i5 - 


-3XA) 


' 8B^C 




A 



(18) 



(19) 



(20) 

Using the above expressions, we can calculate the scalar curvature and the norm of the Ricci 
tensor which are given as, 



Scal = ^(Rc(Fi,Fi)) + (Rc(F2, F2)) + ^{Rc{Fs,Fs)) 

A^X^ + {Bfi-Cu)^ -2AX{Bfi + Cij) 
~ 2ABC 



(21a) 
(21b) 



llRcf = -^(Rc(Fi,Fi))' + -^(Rc(F2,F2))^ + ^(Rc(F3,F3))^ 

_ 3A^X^ - 4A^X%Bfi + Cu) - 4AX{Bfi - Cu)\Bfi + Ci^) + 2A^X\Bfi + Cvf 

(fi/i - Cvf {3B^fi^ + 2BCfiu + 3C^u^) 



(22a) 



4^2 52 C2 



(22b) 
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The higher order flow equations therefore turn out to be, 



^ = -(2Rc(Fi,F0 + a'RcjKTFi)) (23) 
^ = -(2Rc(F2,F2) + a'Rc'(K>2)) (24) 
— = -(2Rc(F3,F3) + a'RcXKTFs)) (25) 
In the next section we analyze the flow equations in five different homogeneous spaces. 

IV. EXAMPLES ON BIANCHI CLASSES 
A. Computation on SU(2) 
1. Flow equations 

The canonical three sphere {S^) is topologically equivalent to the Lie group SU{2), which 
is algebraically represented by 

SUi2) = { :zi,^2eC, \z^f + \z2\' = iy (26) 

All the structure constants are the same in the Milnor frame field. We assume X = fi = 
V = —2. With these values of the structure constants, the 2nd order flow equations turn 
out to be: 

dA _4{B-CY -4A^ _ 1 AA-Bf 
BC "^AB2^ c 

+ (25 + 2A-3C)}2 + -^|^^^^ + (2C + 2A-35)| ] (27) 

dB _ 4{C - AY - 4B^ ^ 1 ^ jC-B)^ 
Itt ~ AC " ^5C^^ A 

+ {2C + 2B- 3A)}2 + J-{iE^^ + (2B + 2A- 3C)Y] (28) 

dC _4{A-B)^-4C^ 1 {A-Cf 

~AB B 

+ {2A + 2C- W)Y + ^{ ^^~^^ + (2C + 2B- 3^)}^] (29) 
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2. Analytical and numerical estimates 



We first try and see if we can analytically estimate the nature of evolution for the scale 
factors A, B, and C. 

(a) Without any loss of generality, we can take the initial values of the scale factors in an 
ordered way i.e. > i?o > Co- Further, from Eqn|27] to |29] we can write the equations for 
the differences (pairwise) of the scale factors as follows. 

ABC 

AA-C) {A' + 2AC + C - B') XA-C)g,{A,B,C) 
ABC A^^ ^^'^ 



dt 




d{A- 


C) 


dt 




d{B- 


C) 


dt 



AB-C) {B' + 2BC + C'-A') ,^, {B-C) GM.B.C) 
ABC A^^ ^^'^ 

where Qi{A, B ,C), Q2{.A, B ,C), Q^^i^A, B ,C) are, respectively 

gi{A,B,C) = {A^ -AA^B + 6A^B^ -AAB'"^ + B^ 

+2A'^C^ + 12ABC^ + 2B^C^ - 8AC^ - 8BC^ + 5C^) (33) 



g2{A,B,C) = {A^ + 2A^B^ -8AB^ + 5B^ -AA^C 

+12AB^C - 8B^C + 6A^C^ + 2B^C^ - 4AC^ + C^) (34) 

gs{A,B,C) = {5A^ -8A^B + 2A^B'^ + B^ -8A^C 

+12A^BC - AB'^C + 2A^C^ + 6B^C^ - ABC^ + C^) (35) 

From the above expressions for the differences, it can be inferred that A{t) > B{t) > C{t) 
holds throughout the evolution, if the initial values satisfy Aq > Bq > Co- Note that each of 
(30)-(32) may be formally solved. For example, if we consider (31), we have A(t) — C(t) = 
(Ao-Co)exp [J)Cit')dt']. 

(b) We now look at the evolution of B{t) for a' = +1. We have, from Eqnl28| 
dB (A^ + C^) ,B^ Pi P2 , (A^ + C^) ,B Pi P2 , , , 
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where Pi = {^^^ + (2C + 2B- 3A)Y and P2 = {^^^^ + (25 + - 30)}^ are strictly 
positive. Choosing Aq = n + 2,Bo = n,Co = n — 2 as the initial values of the scale factors, it 
can be shown that is always negative (from the inequality above, it may be checked 

that the negative term always dominates over the positive term). Further, it is easy to see 
that for n > 4 and a' = (0, +1), all the scale factors are decreasing. 

(c) On the other hand, for a' = —1 the upper bound on B{t) can be written as 

dB B 4M-C)2 Pi P2 

— < -4— + — + 2 — — + 2 — — (37) 

dt - A AC BC^ BA^ ^ ' 

Here we may note a difference with the a' = +1 case. As before, let us choose Aq = 
n + 2, Bq = n,Co = n — 2 to estimate Pi = —4^ + '^^^^q^ — h 2;^2 + 2;^2 which turns 

out to be. Pi = 4(256+n(-64+n(24+n(20-(-3+n)n)))) _ ^^^^ p bouudcd beloW 

(i.e.lim„^oo Pi = —4 and for any finite n. Pi is larger than —4.). Thus, the scale factor P(t) 
can either decrease or increase depending on the choice of initial values. 

(d) It is easy to see from the flow equations that for a' = (0, —1), Aif)^ B{t) and C{t) will 
monotonically decrease when A> [B — C). For typical initial values obeying this condition, 
this decreasing feature is shown in Fig. l(a)^ nd Fig 1(d) where A{t)^B{t) and C{t) are 



obtained by numerically solving the dynamical system. In the graphs, A{t), B{t) and C(t)) 
are shown in thick, dashed and dotted lines, respectively. We follow this convention in our 
plots, throughout this article. 

(e) We now ask whether all the scale factors can be increasing in t. For a' = this is 
not possible. To see this, consider the initial value assumption Aq > Bq > Co- The scale 
factor A{t) will increase for A < (B — C), whereas both B{t) and C{t) will increase for 
A > (B + C) . However, both these inequalities cannot hold simultaneously. Hence all scale 
factors cannot be increasing in t. 

For a' = — 1 we may be tempted to believe in a conclusion similar to the a' = case men- 
tioned above. However, the higher order contribution leads to some unexpected behaviour. 



In Fig, 1(c) we see that the scale factor B(t) initially increases but eventually decreases after 



a certain time. This is due to the the fact that the A > {B + C) criterion does not hold 



in the entire domain of t, a fact apparent from the figure. In contrast, for Fig 1(b) we have 



taken initial values which do not satisfy the aforesaid condition and though, initially, B{t) is 
decreasing we find that after some time, B{t) increases. The reason behind this behaviour 
is — the higher order contribution is much more positive and the net effect makes B{t) 
increase after an initial decreasing phase. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 05 1.0 1.5 2.0 

(c)Ao = 5.5,Bo = 3.5,Co = 1.5,a' = (d) 

-l,Ts = -0.04 ^0 = 7,Bo = 5, Co = 3,a' = 0,Ts = 2.76 

Figure 1: A{t),B{t), C{t) vs t for a' = l,a' = -1 and a' = for SU(2) 

(f) recall that SU{2) admits round Einstein metrics and the Ricci flow converges. This 
feature remains even when we include higher order terms with a' = 1. If we change the 
value of a' (from 1 to 0) a change appears only in the singularity time, though no difference 
appears in the nature of the graph. On the other hand, when a' = —1, we notice an 
expansion of the scale factors - B{t) and A(t) appear to be expanding beyond the minimum 
m Fig. |l(b)| at the same rate, whereas C{t) goes to zero asymptotically. 
Let us now move on to some special cases. 
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3. Special Case: A^B = C 



The flow equations, with this restriction (leading to the so-called Berger sphere metrics 
jit) become, 

dA A^ . A3 

— = -4— - Aa' — 38 
dt ^ ' 

_ = _8 + 4--4.'[-(--l)^ + -(l--)] (39) 

Let us first look at the case, a' = 0. It is clear from Eqn. (38) that A{t) decreases with t and 
has a upper bound, A{t) < Aq — At. However, the same is not true for B{t). If A{t) < 2B{t) 
then B{t) decreases with increasing t, otherwise it increases. B{t) does have a lower bound, 
B{t) > Bq — 8t. Depending on the sign of the R. H. S. of (39) B{t) can have an increasing, 
stationary or decreasing behaviour during its evolution in t. From Eqns. (38), (39), one can 
arrive at a second order equation for B{t) given as: 

,d^B fdB(t)Y dBit) 

There are trivial solutions to this equation: '^'^j*-' = —8 (which corresponds to A{t) = 0) 
and ^^Jj^^ = —4 (which corresponds to A(t) = B{t)). If we define ^ = D{t) then we can 
see the following: (i) At some t if D = Othen^ < (maximum); (ii) If, over a range of 
t, D > then ^ < 0; (iii) If, over a range of t, < then ^ can be < 0, > or = 0. 



We note this behaviour in Fig, 2(d) For example, the maximum occurs at the point where 

D{t) = 0-substitute the values of A{t) and B{t) at the maximum as found in the graph, 

into Eqn. (39) (with a' = 0) and check that it is indeed zero. 

Now we wish to note what happens to the rate of change of {A — B) and {A + B). We 

show that {A — B) decreases faster than {A + B) which implies that as A and B decrease, 

they approach each other. From the above equations for a' = we can write 

d(A-B) -4(A + 2B),^ . ^ A 

^ ' = -{A -B)< - — {A' + 3AB + B^) < -20- (41) 

On the other hand, 

d{A + B) 4 , „ A A 

It is easy to note that A — B decreases faster than A + B. Thus A and B approach each 
other while decreasing. For a' 7^ the analysis can be done following similar logic but we 
prefer to solve the equations numerically and demonstrate our conclusions. 
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(a)Ao = 7, So = 5, a' 1, = 0.914 (b)Ao = 7, = 5, a' -1, -0.602 




(c)Ao = 1, Bo = 0.5, a' = -1, Ts = -0.004 (d)^o = 7, Bo = 5, a' = 0, T, = (-0.76) 
Figure 2: A{t),B{t) vs t for a' = 1, a' = -1 and a' = for SU(2) 

If we plot A(t) (thick line), B{t) (dashed line) treating the evolution equations as a 
dynamical system, we see that both the scale factors converge to the origin (Fig 2(a)"|2(d) ) 
for a' = 0, 1. However, for a' = —1 the evolution of the scale factors depend entirely on the 
initial values of A, B, C- a fact shown in the two figures, Fig.2(b)|2(c 



4-. Phase plots 

We obtain phase plots of the above dynamical system in FigJHl The curves are the 
trajectories of the system for varying initial conditions, and the arrows indicate the direction 
of increasing time. From the phase plots we can deduce the qualitative and quantitative 
behaviour of the flow. We can see that the A = B line is a critical curve demarcating 
the regions in phase space with different behaviour. Also, for a' = all flow trajectories 
converge towards the A = B line. This leads to isotropization of the manifold. However, for 
a' = 1, trajectories from a region in the phase space with A > B converge to the singularity 
B = before isotropization. The behaviour for a' = — 1 is markedly different due to the 
existence of fixed points given by {A,B) = (0, 4), (^, |), (1, 1). In this case, trajectories 
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(a)a' = 1 (b)a' = (c)a' = -1 

Figure 3: 2nd order flow on SU(2)for B = C 



from A > B converge to ( 



16 4 
9 ' 3 



) and those from A < B converge to (0, 1). When A = B, 



we note a convergence to (1,1). Therefore, in this case there is no isotropization, with the 
asymptotic manifold being non-isotropic - a fact which is manifest in the values of the fixed 
points. 

5. Curvature evolution 

Finally, we show the evolution of the scalar curvature. From Eqn. j21b] and Eqn. |22b] 
we can immediately compute the scalar curvature or norm of the Ricci tensor, for locally 
homogeneous three manifolds in the different Bianchi classes. These quantitites depend on 
the scale factors. Thus, if we know the behaviour of the scale factors, under a given flow, we 
can easily obtain the profile of any curvature invariant as a function of the flow parameter. 
However, in most cases, we do not have exact solutions for the scale factors. Therefore, the 
numerical evaluations of the scale factors (for various initial values) are used to obtain the 
curvature invariants. From Eqn J21bl we can see, for SU{2), the scalar curvature is as follows 



Therefore, it can be positive, negative or zero depending on the values of A, B, C . More 
specifically ii A = B = C then F = SA'^ > 0. Similarly if we take A B = C then 
F = A{4B — A), which can be zero, negative or positive for A = AB, A > AB or A < AB 
respectively. 

Fig. 4 shows the profiles for the scalar curvature. For different initial values and three 



Seal 



2 + {B + Cf - 2A{B + C)] _ 2 F{A, B, C) 
ABC ~ ABC 



(43) 
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different ol we have shown the evolution. Both a' = and a' = 1 show a monotonia nature 
of the scalar curvature. For a' = — 1 the scale factors are not always monotonic - the 
behavior depends on the initial conditions. The scalar curvature for a' = — 1 rises initially 
but asymptotically reaches a constant value. The figures in the first row, second column and 
in the third row, second column show the appearance of negative scalar curvature. We note 
here that Seal remains entirely negative or entirely positive during evolution, for a' = 0, 1, 
but may grow from negative to positive value for a' = —1. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.1 0.2 0.3 0.4 




Figure 4: Evolution of Seal for SU{2). 

In summary, several distinguishing features of scale factor evolution and curvature evolu- 
tion arising exclusively due to the presence of the higher order term in the flow, have been 
pointed out in the above analysis and summarised in Table HI 



14 



Table I: Comparison of the evolution of different scale factors and Seal, curvature for SU (2) 





A{t) 


B{t) 


C{t) 


Seal 


a' = 


decreasing 


decreasing 


decreasing 


increases, can be 
-ve. 


a' = 1 


decreasing much 
faster 


decreasing much 
faster 


decreasing much 
faster 


increase /decrease 
(entirely +ve or 
-ve.) 


a' = -1 


depending on 
initial value 


depending on 
initial value 


decreasing 


not monotonic 
(can be +ve, -ve) 



B. Computation on Nil 



1. Flow equations 



v 



Nil is a three dimensional Lie group consisting of all 3 x 3 matrices of the form 

^ 1 X 
1 y 
1^ 

under multiplication, also known as the Heisenberg group. The group action in can be 
written as 

(x, y, z) o (x', y\ z) = {x + x ,y + y , z + z + xy') (44) 
This is a nilpotent group. We can put a left invariant frame field on this Lie group. With the 



inherent metric, fit is a line bundle over the Euclidean plane ¥?. Following Milnor 27(| we 

- — 2 

have the structure constants, A = — 2,yU = = 0. The Rc and Rc tensors can be obtained 
by simple computation. The flow equations in this case will be: 



dA _ 
H " 

dB 
'dt 



4^2 
'EC 

4A 
~C ~ 



-4a' 
20a 



A^ 
BC^ 



(45) 
(46) 



[f] One such left invariant metric is ds^ — dx^ + dy^ + (dz — xdyY 
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dt B 



(47) 



10 - / 

/ 
/ 



5 1/ 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 



(^0, -Bq, Co) 



(a) (b)(Ao,Bo,Co)-(7,5,3),a' 
(7,5,3),a' = l,r, = 0.134 -0.043 



l,Ts 



2 



8 10 



(c) (d)(Ao,Bo,Co) = (3,5,7),a' = l 

(Ao,Bo,Co) = (7, 5, 3), a' = 0,r, = -0.17 

Figure 5: A{t), B{t), C{t) vs t for a' = land a' = for Nil 
2. Numerical estimates 

First we solve the flow equations numerically and look for specific characteristics of the 



flow for different types of initial conditions. In Fig, 5(a) we see that for a' = 1, Aq> Bq> Cq 



and Aq, i?0) Co ^ 1 the flow converges in every direction, contrary to unnormalized Ricci 



flow (Fig. 5(c)) where the 5, C are expanding and diverge diverge while {A{t)) converges 



i.e. a pancake degeneracy exists (the same happens for normalized Ricci flow as mentioned 



m 



12|). Further, for a' = (Fig. 5(b)) and any value of A{t), B{t), C{t), the flow has a past 



singularity. For a' = 1, if Aq < Bq < Co the evolution of the scale factors resemble the 
unnormalized Ricci flow but have no singularity. On the other hand, for a' = —1 (Fig 5(b)) 
and any value of A{t), B{t), C{t), the flow has a past singularity, similar to a' = case. 
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3. Analytical solution 



It is possible to solve the flow equations analytically. Before working out the solutions, 
we note from the flow equations that B = const. C. Let us now choose a new variable 



^ = BC/A. It can be shown easily [3l| that, with this choice and appropriate scaling of 
the coordinates, the nature of the flow can be determined entirely by flnding the evolution 
of ^. In other words. Nil has a one dimensional family of left invariant metrics which is 
parametrised by ^. The flow equations now take the form 

^ = 12 -a' J (48a) 
^ = -Ul + a'l) (48b) 



^ '^l-a'l) (48c) 



where the ^ equation replaces the C equation and A, B are functions of ^. So, effectively, 
there is only one scale factor to worry about, i.e. ^. The ^ equation is readily solved to give 
two solutions- 

C + 3a' \n\^-3a'\ = 12t + k (49a) 
OR 

e = 3a' (49b) 
The equations for A and B can then be solved to obtain A(^) and -B(^)given as, 

(50a) 



A = Aoexp(^-^?j k B = Boexp(^-^?j (50b) 

where (49a) is for the given in (48a) and (49b) corresponds to the solution (48b). The 
evolution of C can also be found using C = In the case a' = 1, and ^ = 3a' we can 
easily see that the scale factors A(t), B(t) and C{t) decrease. When a' = 0, ^ = 12t and 
A{t) ~ t~^, B[t) ~ t^, C{t) ~ t5. For a' = —1, the solution ^ = 3a' is not valid as long as 
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we are dealing with Riemannian manifolds and the evolution is determined from the other 
solution. 

Next we move on to a special case where B = C. 



4. Special Case: A^B = C 

The flow equations, under this assumption are 

dA _ _4A!. _ An' it. 

dt ^" /c-iN 

dt ~ B B3 

It is obvious that the same analytical solutions mentioned above are valid here as long as 
we define ^ = However, we discuss some alternative ways of arriving at the nature of 
scale factor evolution. 

Let us look at the difference between the scale factors given as 

^^^^ = -^^M + B)-Aa\^f{A-W) < -4^^{A-B)-4a'{^^nA-5B) (52) 

It is easily seen that, for a' = (i.e. unnormalized Ricci flow) {A — B) decreases with 
increasing t, but for other values of a' different from zero the scale factors evolve differently. 
The flow equation for B(t) is 

dB f A , A ^r,\ , , 

If ^ < 5 (or < -7) we have 

which explicitly shows that B{t) is decreasing in this region. We then move on to the case 
where ^ > 5. Assume p = The quantity inside the bracket in Eql53]is a polynomial 
in p, f{p) = Ap — 20p^. fij)) represents a parabola (upside down) and its upper bound i.e 
max[f{p)] = I occurs at p = ^ < |. The parabola intersects the abscissa at p = 0, | and 
progressively increases in the negative direction of the ordinate (see Fig|6]). Thus, f{p) is 
positive in the region < p < ^ and, therefore B{t) increases with increasing t. Outside 
this domain of t, B(t) decreases with increasing t. We now show the above-mentioned facts 
numerically. In Fig[|7] we have shown the behaviour of the scale factors for different genres of 



-g2- The first figure (Fig, 7(a) ) has been plotted for such initial values of A{t) and B{t) which 
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0.5 r 




(a)/(^) = 4(^) - 20(^)2 and (^) = i 
Figure 6: /(-^) vs {^) for a' = 1 

correspond to < 5, when both the scale factors are decreasing (see Eqn.[Sl], Fig. |7(b)| ). 



In Fig 7(b) we have shown an interesting turning behavior where B(t) initially decreases but 



eventually increases. The initial condition has been chosen to satisfy > | but during the 
evolution A{t), B(t) decreases so that eventually becomes less than ^ leading to growing 
nature of B{t). On the other hand, for a' = —1,0, the evolution is easy to comprehend 
from the equations (Eqn.[51]). The corresponding evolution of the scale factors is shown, 
respectively in Fig |7(c)| and Fig |7(d)| . 



5. Phase plots 

Let us now turn to the phase portraits(FigJH]). 

The phase portraits in Figj8]are for C = B. From these plots and the equations, it is easy 
to see that A = is a fixed point of the system. For a' = 1, trajectories from one region in 
phase space converge to the singularity 5 = whereas others move towards the fixed point 
A = Q. The critical curve demarcating this behaviour is ^ = BC/A = 3 (plotted in dashed 
line). For a' = and a' = —1 the flow show similar behavior with the trajectories proceed 
towards larger values of B and A = 0. 



6. Curvature evolution 



As before, here also we analyze the evolution of the scalar curvature. We have shown two 
different regimes of initial conditions, namely A^ > Bq > Co and Aq < Bq < Cq in Fig. [9] 
for a' = 1, where we note the increasing or decreasing nature of the scalar curvature (for 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 2 4 6 8 



(a)(Ao, Bo) = (7, 3), a' = 1, T, = 0.03 (b)(Ao, Bo) = (2, 2.45), a' = 1 




2 4 6 8 10 2 4 6 8 10 

(c)(^o,5o) = (7,1), a' = -l,r, = -0.0003 (d)(Ao,Bo) = (5, 3), a' 0,r3 = -0.15 
Figure 7: A(t),5(t) vs t for a' = 1, a' = -1 and a' = 




(a)a' = 1 (b)a' = {c)a' = -1 

Figure 8: 2nd order flow on Nil for B = C 

different initial conditions). For a' = and a' = —1 the scalar curvature monotonically 
converges to zero asymptotically after beginning from a negative value. 

We conclude this section by providing a comparison of the evolution of different scale 
factors and Seal, curvature for Nil manifolds in tabldlTl 
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(a)Evolution of scalar curvature for various a', and Aq, Bq, Co 
Figure 9: Evolution of scalar curvature for Nil. 



Table II: Comparisons of the evolution of different scale factors and Seal, curvature for Nil 





A{t) 


m c{t) 


Seal 








asymp. flat 


a' = 


decreasing 


increasing increasing 










starting from -ve. 








increase/decrease 








entirely -ve, 


a' = 1 


decreasing 


depends on i.c depends on i.c 










past / future 








asymp. flat 




decreasing, 




asymp. flat 


a' = -1 




increasing increasing 






asymp. flat 




starting from -ve. 


C. Computation on Sol 






1. Flow equations 






Sol is a solvable j 


j;roup which has 


minimum symmetry among all the 


eight geometries. 


The group action in 


can be written as 





(x, y, z) o (x', y', z') = {x + e ^x', y + e^'y', z + z') 



(55) 
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We can put left invariant vector fields^ on the manifold, for which the structure constants 
in the Milnor frame will be A = —2, /i = 0, = +2. Sol really has a two parameter family of 
metrics upto diffeomeorphism (obtained by scaling and redefinition of the coordinates, see 



3l|). Using the values of the structure constants, we can find the components of Rc and 



— 2 

Rc tensors. The 2nd order flow equations reduce to. 



11 = BC A^- ^^^^ 



dB _ {A + Cf {A + Cfi^A^ - QAC + 5C^) 

dt AC A^BC^ ^ ' 

dC ^ jC' - A^) _ , {A + Cn5A'-2AC + C') 

dt AB A^B^C ^ ' 



2. Numerical and analytical estimates 

(a) The behaviour of the scale factors for the above higher order flow are different from those 



in unnormalized Ricci flow. From Fig 10(d) (unnormalized Ricci flow) we note that the scale 
factors do not have a future singularity. In fact, we see the appearance of a cigar degeneracy. 
The above behaviour can be predicted qualitatively from the equations themselves. We note 
that A and C can be interchanged in the equations. So instead of three equations it is 
enough to examine two of them. Further, without loss of generality, we may assume A > C. 
For a' = 1 the scale factors converge (see Fig 10 (a)). If a' = and A > C we find that B{t), 
C(t) increase and A{t) decreases (more detail for a' = can be found in [l|]). If a' = — 1 the 
behaviour of the scale factors bear a resemblance with unnormalized Ricci flow, though the 



singularity time (in the past) changes due to the higher order term ( Fig, 10(c) illustrates 
one such example for a particular set of initial values). 

(b) We now consider the case a' = 1 and analyze the evolution of the scale factor A{t) a 
bit further. It can be shown that Eqnj56]can be written as 

dA .AA(^ 2A(^ CV AV %A(^ ^ 



which shows that Aif) is decreasing in forward time. This is also true for Cit) where 
the second term in Eqj58] dominates over the first term and, therefore, the net effect is a 



[^] it is easy to check that one such left invariant metric will be ds^ = e^^dx^ + e ^^dy 
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decreasing C{t). In the same way, we can argue that B{t) is increasing (decreasing) for 
A < B [A > B). Since A{t) and C{t) are both decreasing, they may approach each other. 
To show this, consider the difference A — C, which satisfies 



d{A-c) M-c) 



dt 



-4 
< -4 



ABC 
{A + Cf 

ABC 



{A + Cf -4a 



,{A + CYiA-C) 



A^B^C^ 



{A^ - 6AC + C^ 



C - 16a 



>{A + Cf 
B^C 



(60a) 
(60b) 



and demonstrates that the rate of change of A — C is always negative. This feature is also 
clearly visible in Fig. 10(b). If a' = 0, we may note that ^ < -AA or A{t) < {Aq - 4t). 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.0 0.2 



(a)(Ao,Bo,Co) = (7,5,3),a' = l,r, 
0.078 



(b)(Ao,Bo,Co) = (5,7,3),a' = l 



40^ 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

{c){Ao,Bo, Co) = (7, 5, 3), a' = -1, = (d)(Ao, Bq, Cq) = (7, 5, 3), a' = 0, T, = 
-0.052 -0.154 

Figure 10: A{t), B{t), C{t) vs t for a' = 1, a' = -1 and a' = in Sol manifold 



Next we move on to a special case where A = C, which is analytically solvable. 
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3. Special Case : A = C 



With this assumption, the flow equations simphfy, and we have 

dA , A 

f^^ = 8-32a'i 

dt B 



(61) 



(62) 



When a' = 0, A(t) is a constant and B{t) increases hnearly. If a' = 1, A(t) decreases 
but B{t) may increase/decrease depending on initial conditions (as shown in Fig 11 (a), (b)). 
For a' = — 1, A and B both increase in forward time (Fig. 11(c)). 

The flow equations can be easily solved to obtain explicit solutions as - 




(a)(^o,So) = (3,7),a' = l 



10 0.00 0.05 0.10 0.15 0.20 0.25 0.30 



(b)(Ao,Bo) = (7, 3), a' = l,r, = 0.318 



20 - 



ic){Ao,Bo,Co) = (3, 5, 7), a' = -l,r, = 
-0.219 

Figure 11: A{t), B{t) vs t for a' = I and a' = -1 
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(a)a' = 1 (b)a' = (c)a' = -1 

Figure 12: 2nd order flow on Sol for A = C 



( 4a' \ 

A\\ - — \ = constant & B + ia' In \B - 4a'| = 8t + k 




(63a) 



OR 



A = Aoexpi -t\ & B = Aa' 



(63b) 



Note that the second solution is not valid as long as we are concerned with Riemannian 
manifolds. 

4- Phase plots 

The phase space is plotted in FigJT2l For a' = 1 flow trajectories from B < 4a' converge 
to the singularity B = whereas others tend towards the flxed point A = 0. B = 4a' is 
the critical curve in this case. For a' = flows we have A = Aq and B = 8t + Bq, with the 
trajectories being straight lines flowing towards ever increasing values of B. In the a' = — 1 
case the trajectories expand outwards from A = B = to larger and larger values. 

5. Curvature evolution 

Lastly we examine the evolution of scalar curvature. For a' = and a' = — 1 the behavior 
is similar to that of Nil. The only difference appears in the a' = 1 case where we see that 
the curvature increases for Cq < Aq < Bq while it decreases in the other regime. 
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(a)Evolution of scalar curvature for various a', Aq, Bq, Cq 
Figure 13: Evolution of scalar curvature for Sol. 



D. Computation on Isom(M2) 



1. Flow equations 



This case also involves a kind of solvable geometry. It is not mentioned in Thurston's 
eight geometries because it is not maximally symmetric. According to Milnor's prescription 







somfM^) also, like Sol 



31|). As done earlier, we 



we have the structure constants as follows A = —2,fi = —2,iy 
has a two parameter family of metrics upto diffeomeorphism (see 
can compute the various curvature quantitites and end up with the 2nd order flow equations 
given as 



dA _ 
Itt ~ 
dB 

It 



-A{A' - B') 

5C 
-A{B^ - A^) 

~AC 
A{A-Bf 

AB 



- 4a 

- 4a 
4a 



,{A- Bf{A^ + 2AB + ^B^ 



AB^C^ 

,{A- BfihA^ + 2AB + B"^) 
A^BC^ 

,{A - BfihA^ + QAB + hB"^) 
A^B'^C 



(64) 
(65) 
(66) 



2. Analytical estimates 



(a) The Isom(M2) class does admit Einstein metrics (flat metrics for this case) and converges 



under unnormalized Ricci fiow(Fig, 16(a) ). It is easily observable that A = B is the fixed 
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point of the flow irrespective of the presence of the higher order term. For a' = 1 each of 
the scale factors attain constant value asymptotically and they are seen to decrease initially, 
(b) Let us check the evolution of each of the scale factors for a' = 1. Using the fact that 
{A^ + 2AB + 5B^) > {A + and EqnElwe can write the evolution of A{t) as 

'^A < M-Bf _ M + mA-Bf 

dt - BC AB^C^ ^ ' 

So the scale factor A{t) always decreases. Using the same kind of argument as earlier 
(5^2 + 2AB + 5^) > {A + BY we can rewrite EqJSJ] as 

dB AA + Bf _ {A + Bf{A-Bf ( {A-Bf \ {A-Bf 

dt - AC A^BC^ - V ^BC J AC ^ ^ 

It can be easily shown that ^1 — ^"^^-^ (say Fi) can be negative as well as positive. To 
illustrate this we assume A = n + 2,B = n,C = n — 2 and calculate Fi which turns out 
to be 1 — („_2)ra(2+n) ' ^rom the plot of Fi (Fig JT4|) . we note that depending on the value of 
n, Fi can be positive as well as negative. Obviously this holds for B{t) as well. Next, we 
estimate the behavior of From EqnEHlit is easy to anticipate that (SA^ + QAB + 55^) 
will be the deciding factor. Using the bound — (5A^ + QAB + 55^) < —{A — B)^ we can 
recast the Eqnj66] as 



dt - 'B ' \A AC ) 
The behaviour for may be found by choosing A, B, C and calculating a quantity similar 
to the Fx mentioned earlier. The evolution of Cit) is qualitatively similar to that of -B(t). 

3. Numerical estimates 

In Fig. [15] we demonstrate our conclusions for certain specific inital values by numericaly 
solving the dynamical system. The cases with a' = 1 and different sets of initial values are 



shown in Fig. 15(a)-(c). For a' = —1 the flow develops a past singularity (Fig 15(d)). The 



nature of the evolution of scale factors for a' = is similar to the case a' = — 1 except for 



singularity time. These features appear in Fig. 16(a) and Fig, 16(b) 



J^.. Special case: B = rjA 

Let us now choose B = rjA with < ?7 < 1. The flow equations become - 
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Figure 14: Fi[n) vs n 
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-0.128 

Figure 15: C(t) vs t for a' = 1 and a' = -1 in Isom(M2) 
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{a)(Ao,Bo,Co) = (7, 5, 3), a' = 0,r^ = -0.3 (b)coniparison, a' = -1, a' = O(thick) 



Figure 16: A{t) , B (t) , C (t) vs t for a' = in Isom(R2) 



dA _ 4^ (^^ - 1) _ 4 / MV - l)2(5r/2 + 2r/ + 1) 
Chi " cv 



(70a) 



(70b) 



dC_^ {il~lf (^-1)^(5//^ + 6// + 5) 



(70c) 



For a' = we can exactly solve the 7] — C system and obtain a relation between r]{= ^) 
and C which is given as 



m = B{t) 



(71) 



(^k + ^k^ - 4C{t) 

We substitute the above relation back in EqnJTOb and arrive at an exact solution for r]{t) as 



2v^ +l^(v^-l) 



32t 



(72) 



(1 + r/) (v^+1) A; 
We have analyzed numerically the r] — C system in FigJTTl We note that r/ = 1 is a fixed 
line and attracts fiow trajectories in all cases. For a' = 0, — 1 all trajectories in phase space 
end up on ?7 = 1. However, for a' = 1 some fiow to 77 = 1 whereas others fiow towards the 
singularity at 77 = = C. 
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(a)a' = 1 



(b)a' = 



(c)a' = -1 



Figure 17: 2nd order flow on Isom(IR^) 
5. Curvature evolution 

Before we end our discussion, it is worthwhile to mention the evolution of scalar cur- 
vature for this case. Like the Sol and Nil manifolds, the scalar curvature, for a' = 0, — 1, 
asymptotically reaches zero but by starting from a negative value. The case a' = 1 has 
different behavior depending on the initial conditions, as seen in the Fig. [TB]. For some 
particular initial values, the scalar curvature starts from a negative value, decreases first, 
then asymptotically reaches zero. For the same a' = +1, but different initial conditions, we 
obtain negative scalar curvature as before, but it grows to larger negative values till the flow 
exists. 

E. Computation on SL(2,M) 
1. Flow equations 

This is a three dimensional Lie group consisting of all 2 x 2 matrices with unit determi- 
nant and its universal cover is denoted by SL(2,]R). It can be shown that SL(2,M) is not 
isometric to x M which is also a locally homogeneous space. We can put a left invariant 
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(a)Evolution of scalar curvature for various a', Aq, Bq, Cq 



Figure 18: Evolution of Scalar curvature for Isom(M2). 

metric| (or right invariant metric but not bi-invariant) on this manifold. Using Milnor's 

prescription (A = —2, /i = —2, u = +2) we find the various curvature-related quantities. 

^ — 2 

The evaluated components of Rc and Rc eventually lead to the following flow equations: 



dA _ ^A'^-{B + Cy ^^,2{A^ + 2A^{B + Cf -8A{B - C){B + C) 



dt BC B^C^A 

,2{B + C)2(5fi2 _ QBC + 5C2) 



-2a 



B^C^A 



dB_ _ B^ - {A + Cf _ , 2(5A^ + AAC{B + Cf + ^^(-8^ + AC)) 

'dt ~ AC " AKPB 

,2(2^2(^2 _ 4BC - C^) + {B + C)\B^ - 2BC + 5C^)) 



-2d 



A^C^B 



dC 2C'^ - 2{A- Bf ^ ,2{hA^ -AA^[^B-2C)-AAB{B^Cf 
^ -2 : 2a 



dt AB A^B'^C 

,2{-2A^{B^ + ABC - C^) + (5 + C)^{5B^ - 2BC + C^)) 

A^B^C 



(73) 



(74) 



(75) 



[I] one such set of left invariant forms are 

TT^ = + (1 + x'^)'^dy + (x — y — x'^y)dz 
< TT^ = 2xdy + (1 - 2xy)dz 
^ TT^ = dx + {x^ — l)dy + {x + y — x'^y)dz 
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2. Analytical and numerical estimates 



(a) The abovementioned equations are very hard to solve, analytically, when a' 7^ 0. For 
a' = we can see from Eqn|7n]and EqnJTHthat the flow is symmetric in A and B. Assuming 
^0 ^ -So, without any loss of generality, we can show that Aii) > B(t) throughout the 
a' = flow. This may be inferred (for A > B > C) from the evolution of the difference of 
{A — B) given as 

'^-^^^ = -^{A-B){A + B-C){A + B + C)<Q (76) 

The evolution of C{t) is straightforward. We note that there exists a lower bound, in 
the following sense. 



'^ = ^UA-Bf~C') 
dt AB^^ ' ' 



4( 



A B 



2) > -4 



(77) 



B A) AB 

EqnJTTl shows that C(t)(= Cq — 4t) is monotonically decreasing -a point of difference from 
what we find for normalised Ricci flow. Similarly, it is not difficult to show that -B(t) is also 
monotonically increasing, if we assume > -Bo > Co- flowever, the nature of the evolution 
of A{^^ shows an increase though it is not monotonic. This can be justified as follows. We 
have 



dA 

H 



BC 



{A + B + C){B + C ~A) 



(78) 



Thus, A{t) will increase monotonically provided A < {B + C), otherwise it will decrease 
initially and then increase. We further note that as A{t) and B{t) increase they approach 
each other. This feature follows from EqnJ76l assuming Aq > Bq > Cq. We can write 



d{A - B) 
Jt 



< 



{A- B){A + B-C){A + B + C) 



ABC 

^ {A- B){A + B -C) 



AB 



< -12- 



C 
AB 



(A-B) 



All the above stated features for a' = are shown in Fig, 19 (e 



(79) 
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(b) For a' 7^ it is difficult to understand even the qualitative nature of the evolution of 
scale factors, so we tried numerical solutions for some particular initial values. Before we 
start discussing the numerical results, let us recall that the SL(2,M) class does not contain 
Einstein metrics. The unnormalized Ricci flow does not converge- instead it evolves towards 
the pancake degeneracy. But the inclusion of higher order terms seems to oppose this. For 
a' = 1 the scale factors converge for different initial values -a fact depicted in Figs. (19)(a-c). 
On the other hand, the a' = —1 case resembles the unnormalized Ricci flow where two of 
the expanding scale factors approach each other, though the past singularity time is larger 



(Fig, 19(d) ). The generic behavior of scale factors for a' = are retained for a' = —1, except 
for the evolution of C(t)-a fact depicted in the Figs. (19)(e-f). 

Finally, we move on the special case where A = B ^ C where we have an exact relation 
between the scale factors for a' = 0. 



3. Special case: A = B ^ C 

The flow equations in this case are: 

(80b) 

For a' = we can exactly solve the system and arrive at a relation between A and C given as 

^ ^ A{t)' h expi-Ajif) 
§1 + fci exp(-A(t)2) 

We can also check that for a' = the scale factors diverge from each other. This may be 
noted from the rate of evolution of the difference of scale factors, as given below. 

d(A-C) C f C\ 

-°'- l + - >8 (82) 



dt A \ A^ 

For a' = — 1 we can write down the expressions for the differences of the scale factors from 
EqnJHUl as follows, 

d{A-C) ^ 4(5C2 + A2(8 + C) + AC(12 + C)) 4 /CV 

+ T T y^-^) 



dt A3 A \A 

I 4(A + C)(5C + A(C + 7)) ' ^g^^ 
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(a)(Ao,-Bo,Co) = (9,7,5),a' = 1 (b) 

(^0, Bo, Co) = (3, 5, 7), a' = 1, = 0.033 





0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.2 0.4 0.6 0.8 1.0 



(c)(v4o, Bo, Co) = (7, 5, 3), a' = 1, = 0.6 (d) 

(Ao,Bo,Co) = (7,5,3),a' = -l,T, = -0.07 




0.0 0.2 



1.0 0.0 0.2 0.4 0.6 0.8 



(e) (f)a' = -1 and a' = (thick) 

(^0, Bo, Co) = (7, 5, 3), a' = 0, = -0.24 

Figure 19: A{t), B{t), C{t) vs t for a' = I, a' = -land a' = for SU2^) 
where we maintain A > C . Here the scale factors are diverging. This fact has been nu- 



merically verified in Fig, 20(c) We do not show the evolution of the scale factors A{t) and 
C{t) for a' = because their nature are similar to the a' = —1 case. Lastly, if a' = 1, the 
evolution of the scale factor C{t) can be found from Eqn JSObl However, the evolution of A{t) 



will be different for different initial conditions — features which are shown in Fig, 20(a) md 
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inFig, 20(b) Let us now turn to the phase potraits. 



4. Phase plots 

As before, we plot the above-mentioned reduced system in Figj2TJ For a' = 0,-1 tra- 
jectories flow to larger values of A. But in the case of a' = 1, the presence of a fixed point 
{A, C) = (4a', 0) makes the flow more interesting. In this case, some trajectories from a 
certain region in the phase space go towards larger A as before, but others converge to the 
singularity at A = 0. There also exists a critical curve, which flows into the fixed point, and 
demarcates the regions with different asymptotics. 
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(a)(ylo,Co)-(7,5),a' = l 



6 0.00 0.02 0.04 0.06 O.OS 0.10 



(b)(Ao,Co) = (3.5, 2), a' = l,r, = 0.11 




(c)(Ao,Co) - (7, 5), a' = -1,T, = -0.13 
Figure 20: A{t), C{t) vs t for a' = 1, and a' 



5. Curvature evolution 



We end by plotting the scalar curvature evolution. This is not very different from the 
earlier cases (Sol and Nil), for a' = and a' = — 1 where the curvature asymptotically goes 
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(a)a' = 1 (b)a' = (c)a' = -1 



Figure 21: 2nti order flow on SL(2,M) ior A = B 

to zero starting from a negative value. However, when a' = +1 we find some interesting 
features which can be seen in Fig. [22] . For a particular initial condition, the scalar curvature 
increases from a negative value towards zero, stays there for a while and then drops to lower 
values once again. 




02468 10 02468 10 



(a)Evolution of scalar curvature for various a', Ao,Bo,Co 
Figure 22: Evolution of scalar curvature for SL(2,M). 
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V. CONCLUSIONS 



Our overall aim in this work has been to study in detail the consequences of second 
order (in Riemann curvature) geometric flows on three dimensional homogeneous spaces, 
using analytical, semi- analytical and numerical methods. Through the analysis carried out, 
we believe we have been able to obtain quite a few characteristics which seem to arise for 
second order geometric flows on three dimensional homogeneous (locally) geometries. Here, 
we briefly summarize our results and mention a few possibilities for the future. 

From our results, we can say that for manifolds which do not contain Einstein metrics the 
flow characteristics show new features essentially caused by the inclusion of the higher order 
term. On the other hand, the class of group manifolds which admit metrics of Einstein class 
(whether it is flat or round), results for second order flows appear to be refinements over 
known results for Ricci flows. However, generic new characteristics do arise with varying 
sign of OL -in particular, a negative a'- and these have been noted in our work. 

In several restricted cases (i.e. where two of the scale factors are related), we are able 
to solve the flow equations exactly. We have obtained analytical solutions in such restricted 
cases for Nil^ Sol, Isom(M2)^ SL(2,M) manifolds. The exact solutions are instructive because 
they help in obtaining analytical expressions for fixed points (curves) as well as in under- 
standing the evolution of the scale factor. In addition, we also use them for checking our 
numerics. 

A generic observation is the fact that the singularity time changes due to the inclusion of 
higher orders in the flow equations. This pattern is noticeable throughout in our numerical 
work. 

The results for a' < are, in quite a few cases, strikingly different from those for a' = 
or a' > 0. Even the evolution of the scalar curvature exhibits a different behaviour in many 
of the cases studied. 

For SU{2) all the scale factors do not converge for a' = —1 which is exactly opposite to 
the characteristics for a' = 1 and a' = 0. Here, if a' = 0, 1, the scalar curvature increases 
but when a' = — 1 it increases first and then decreases. The appearance of negative scalar 
curvature in SU{2) is also noted. In the case for Sol manifold the behavior of the scale 
factors depend on different initial conditions for a' = ±1 — so does the evolution of the 
scalar curvature. 
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In the case of Isom{M?) all the scale factors converge for a' = 1. When a' = —1 one 
scale factor decreases while the other two increase. In the last case, SL(2,R), we find that 
for a' = 1 all the scale factors may initially increase but they converge towards a singularity. 
However, for a' = —1 all the scale factors diverge. 

In all cases we have obtained the phase portraits, which we feel, helps in visualising the 
flow features as well as the fixed points. We also provide a summary of all our results in 
several tables in the sections as well as, in the end. 

Among possible future directions, we mention a few below. 



13| for such higher order 



• It would be interesting to pursue the approach presented in 
flows. 

• A more complete analysis of the stability and classification of fixed points can be worked 
out for such second order flows. 

• Given the fact that such three manifolds do arise in various physically relevant contexts, 
it will be nice to know whether our results on higher order flows can help us understand 
such scenarios in any meaningful way. 

• The relevance of our results in the context of renormalisation group flows of the bosonic 
nonlinear a-model deserve some attention. 

• Since homogeneous four manifolds are already classified and studied with reference to 
Ricci flows 32| it would be worthwhile extending our results to four manifolds. 

We hope to address some of these issues in future articles. 
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Table III: Comparisons of scale factors for different cases 



Manifold Bianchi Type 


q' = 


a' = 1 


a' = -1 


SU(2) 


IX 


all scale factors converge to finite 
time singularity 


different convergence rates, 
singularity times 


does not always converge to a 
singularity 


Nil 


II 


pancake degeneracy 


depends on initial data 


pancake degeneracy 


Sol 


IV- 1 


cigar degeneracy 


depends on initial data 


depends on initial data 


Isom(R)2 


VIIo 


one decreases and other two increase, 
converges to flat 


all decrease and converge to flat 


one decreases and other two increase, 
converges to flat 


SL(2JR) 


VIII 


one decreases and two increase 


increases or decrease 


all increase and diverging 



Table IV: Comparisons of scalar curvature between different cases 



Manifold 



a' = 



a' = 1 



a 



SU(2) 



Seal, curvature 



Seal, curvature 



Seal, curvature 
increases can be -ve. 



increases or decreases increases/decreases, 



(entirely -ve. or +ve.) 



can be -ve. 



Nil 



asymptotically goes to 
flat from -ve. value 



increases 
(^0 < ^0 < Co) or 
decreases 
(Ao >Bo> Co) 
depending on initial 
value 



asymptotically goes to 
flat from -ve. value 



Sol 



asymptotically goes to 
flat from -ve. value 



increases 
(Co <Ao< Bo) or 
decreases 
{Ao >Bo> Co) 
depending on initial 
value 



asymptotically goes to 
flat from -ve. value 







starts with -ve. value, 




starts with -ve. 


value, 


may develop 


starts with -ve. value, 


Isom(IR)^ asymptotically 


goes to 


singularity. 


asymptotically flat 


flat 




asymptotically flat 








depending on initial 




^-^^ asymptotically 


goes to 


value, curvature may 


asymptotically goes to 


SL(2,E) 








flat from -ve. 


value 


decrease (tend to zero) 
or increase 


flat from -ve. value 
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